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Abstract. We prove several claims made by Kontsevich about the orbifold 
Euler characteristic of the three types of graph homology introduced by him. 
For this purpose, first we develop a simplified version of the Feynman diagram 
method, which requires integrals in one variable only, to obtain the generat- 
ing functions as asymptotic expansions of certain Gaussian integrals. Finally, 
following Penner, we relate these integrals to the gamma function in order to 
compute the individual coefficients of the generating functions. 



1. Introduction 

In two papers ([10], [11]) published in 1993 and 1994, Kontsevich constructed 
a family of objects called graph complexes. These are chain complexes of vector 
spaces, where each vector space is spanned by a set of finite graphs equipped with an 
orientation. He computed the orbifold Euler characteristic of these chain complexes, 
and stated the answers; however, he did not give any proofs except for a one sentence 
reference to Feynman diagrams. 

In Sections 5, 6 and 7 we recall the definitions of the graph complexes and their 
orbifold Euler characteristic, and prove most of Kontsevich's claims about them. 
But first, we introduce the two technical tools used in the proof. In Section 3, 
we develop a simplified version of the Feynman diagram apparatus (also called 
the method of Gaussian integrals). This method, well known to physicists but 
mostly unknown to mathematicians, can be used to construct generating functions 
for various problems involving the counting of graphs. A good exposition of the 
general method can be found in [4]. Our version is based on [15]. It is hmited in 
that it cannot take into account the number of boundary components of ribbon 
graphs (referred to as the "topological expansion" in the physics literature), but 
it is much simpler as it requires integration over M only, instead of over a space 
of matrices. In Section 4, we recall some basic facts about asymptotic expansions, 
which are used to extract information about the individual terms of the generating 
functions obtained by the Feynman diagram method. 

The emphasis is on proving the claims made by Kontsevich about the orbifold 
Euler characteristic of the commutative, associative and Lie graph complexes. How- 
ever, our simplified Feynman diagram method can be applied to graph complexes 
based on other cyclic operads (see [5] or [13]), as well as to other problems involving 
the counting of graphs. 
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3. Enumeration of graphs 

As our first example, we consider the following counting problem. Let T2e be 
the set of graphs with e edges such that each vertex has valence at least 3. (The 
valence of a vertex is the number of edges incident to it.) For each G £ T2e, let 
Aut(G) be the set of automorphisms of G. The goal is to evaluate the sum 

The first observation is that we are, in fact, counting the number of labeled 
graphs. We will say that a graph with e edges is labeled if the 2e half-edges are 
numbered from 1 through 2e (see Figure 1). Fix a graph G, and consider all its 
labelings, i.e., the (isomorphism classes of) labeled graphs with G as the underlying 
unlabeled graph. The symmetric group S2e acts on these labelings, and the size of 
the stabilizer of any given labeling is equal to the number of automorphisms of G. 
Thus if T|g denotes the set of labeled graphs with e edges such that each vertex is 
at least trivalent, then 



(3.2) 



1 _ ^ I {labelings of G} I _ ^ 1 



^ |Aut(G)| ^ (2e)! ^ (2e)! 



i.e., the sum to be computed is the number of labeled graphs jTjgl divided by (2e)! 




Figure 1. Labeled graph = chord diagram + partition 

Second, notice that we can break up a labeled graph into two separate objects: 
a pairing of the numbers 1, 2, . . . , 2e, which we will call a chord diagram, describes 
the edges; and a partition of these 2e numbers into groups of three or more tells 
us what the vertices are. Conversely, any pair of a chord diagram and such a 
partition corresponds to a labeled graph. Therefore, if C/i„ denotes the number of 
chord diagrams on n numbers and Par^ the number of set partitions of n distinct 
objects into subsets of size three or more, then 

(3.3) \Tl\ = Ch2eParl+. 

Each of the quantities on the right hand side are fairly easy to compute. In the 
case of chord diagrams, this could be done directly, but in both cases the simplest 
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solution is to find their exponential generating functions Ch{x) = Yll^=o Chnx"/nl 
and Par^+{x) = ^^=0 Parl+x^'/nl. 

Namely, let Q denote some class of objects; following [9], wc will call it a species. 
Let Q{x) = ^'^^qQux" /nl be the corresponding exponential generating function, 
where Qn is the number of objects of type Q on n numbers, vertices etc. Consider 
the exponential species exp{Q), the objects of which are the disjoint unions of any 
number of objects of type Q. Then the exponential generating function {exp{Q)){x) 
is equal to the formal exponential 

(3.4) exp(Q(a;)) := 1 + Q{x) + ^{Q{x))' + ■■■ + ^{Qi^))" + ■■■ 

of the power scries Q{x), provided that any object of type exp{Q) can be uniquely 
decomposed into its sub-objects of type Q. In particular, we must have Qq = 0, so 
the sum above is finite in each coefficient. This is an application of Joyal's theory 
of combinatorial species; details can be found in [3] or [12]. 

Now, let E"^ be the species of two-element sets, i.e., _E| = 1 and = for all 
n 7^ 2. The exponential generating function is E'^{x) = x^ /2. On the other hand, 
we have Ch = exp{E^), and thus 

(3.5) C/i(a;) =e"'/2^ 

Similarly, if denotes the species of sets of size greater than or equal to three, 
then we have 

(3.6) E'+{x) = ^ + ^ + .--=e^-l-x-x'/2. 
Since Par^+ = exp(£^^+), we immediately obtain 

(3.7) Par^+{x) = e'^'^-^-^-^'/^. 

Finally, we would like to find the exponential generating function of |T^|, given 
those of Ch and Par^^ , using equation (3.3). This is called the Cartesian product 
of the two species, denoted by |T^| = Ch x Par^^ . There does not seem to be 
any general method for computing the generating function of a Cartesian product; 
however, in this particular case, when one of the factors is Ch, the following lemma 
solves the problem. Readers familiar with the the theory of Feynman diagrams 
will recognize it as a one-dimensional version of the calculation underlying Wick's 
Lemma. 



Lemma. C/i(x) = -j= / e'^^ '"^ dy . 
\/2tt J -oo 



Proof. Straightforward. □ 

Another way of stating the lemma is 

(3.8) Chn = \ e-^'/2yn fQj. ^_ 

\/27r J_oo 

Note that this formula holds for odd values of n, as well: in this case, both sides 
are zero. If we could exchange the order of the summation and the integration, we 
would have 



(3.9) \T\x) = ^\ e-^'/V^'-i-^^-^'^'/^rfy. 

V 27r J-oo 

Although this formula needs explanation (the series |T"^|(x) does not converge: 
its radius of convergence is zero), it is indeed valid in a limited sense, as we shall 
see below. 
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4. Asymptotic expansions 

Power series with zero radius of convergence can be treated analytically by con- 
sidering them as the asymptotic expansion of some function. We briefly recall the 
definition and basic properties of asymptotic expansions; details can be found in 

[6], for example. 

Definition. Let / be a function defined on a subset of K including a (possibly 
one-sided) neighborhood of 0, with values in M or C. We say that the formal power 
series X^^q '^n^" asymptotic expansion of / around a; = and we write 



(4.1) f{x) w ^2 ^n^" as a; ^ 

n=0 

if the following condition is satisfied for each N > 1: 

N-l 

(4.2) f{x) = J2 o-nx"" + 0{x^) as x 



n=0 



This means that for each fixed A'', the quotient {f{x) — Y^'^=Q anX")/x^ is bounded 
in some neighborhood of 0. 

Note that the neighborhoods wo take in (4.2) may shrink down to zero as TV ^ oo, 
in which case our definition does not say anything about the tail of the series 
X^^o (^nX"" for any fixed non-zero x. It is possible that f{x) « J2'^=o ^-nX^, but 
the scries docs not converge for any non-zero value of x. Even when the series 
does converge, its sum is not necessarily equal to f{x); for example, e~^/^ « + 
Ox + Ox^ + • • • as X — > 0+. Also note that the function / does not need to be 
diffcrcntiablc, or even continuous. 

On the other hand, if / is a C°° function, and its Taylor series has a positive 
radius of convergence, then the Taylor series is an asymptotic expansion of /. A 
function may not have an asymptotic expansion, but if it docs, it is unique: one 
can show by induction on (4.2) that the two series are equal term by term. All 
the usual operations on formal power series carry over to asymptotic expansions: 
given the asymptotic expansions of / and g, the functions / + g, fg, eJ etc. all 
have asymptotic expansions, which are equal to the sum, product, exponential etc. 
of those of / and g, as formal power series. 

If we consider it as an asymptotic expansion, any power scries can be exchanged 
with the integral in the Lemma of the previous section. More precisely, the following 
is true. 

Proposition. Suppose that ft,: M ^ C has the asymptotic expansion around zero 
h{x) X^^o '^nx'^, and the integral f{x) = /'^h{xy) dy converges for 

all X in some neighborhood of zero. Then f has the following asymptotic expansion 
around zero: 

oo / -. poo \ 

/(x)«V^/ e-y"^^y"dy)anx" as x ^ 
„_n V v27r J-oo J 

« Chnttnx'^ as X — > 0. 

n=0 
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Proof. Since h{x) = J2n=o + 0{x^), we have 

(4.4) -= / e-y'/'h{xy)dy 

(4.5) = / e-y-/'' + '^(^'"y") dy 

(4.6) = ^ / V a„x"y" dy + ^ e-y' ^'Oix" y") dy 
„=oVV27ry_oo / 



(4.7) 



because e~y^ I'^y^ dy is finite for any fixed value of A''. □ 
Or in other words: 

Theorem 1. Let P he any species. If there is a function /i: R — > C such that P{x) 
is the asymptotic expansion of h around zero, then 

1 1"^ 2 

(4.8) -= I e-y^l'^h{xy) dy w {Ch x P){x) as a; ^ 
V27r J-oo 

as long as the the integral converges for all x in some neighborhood of zero. 

We will also need the following more general version, which can be proved in the 
same way. 

Corollary. Let h: WxM. ^ C be defined on U xM., where U is some neighborhood 
of zero, and suppose there is a sequence Pn of polynomials such that 

N-l 

(4.9) h{x, y)=J2 Pn{y)x'' + 0(x^y^(^)) asx^Q 

n=0 

holds for each N > 1, where M{N) is an integer depending only on N. Suppose 
further that the integral f{x) = /^h{x, y) dy converges for all x G U. 

Then f has an asymptotic expansion around zero, which can be obtained from (4-9) 
by replacing each occurrence of y^ by the number Chk- 

5. The orbifold Euler characteristic 

Our goal is the compute the orbifold Euler characteristic of certain graph com- 
plexes. A graph complex is a chain complex 

• • • > Ofe+i > Ofe > > • • • > Oi > Oq — <■ U 

where each Cfe is a vector space (in this paper, always over M) spanned by some set 
of graphs. In our examples, graphs in Cfc have k vertices. Given a graph complex 

C* of finite length where each Ck is spanned by a basis corresponding to the finite 
list of graphs Gfei, Gk2, ■ ■ Gknk^ its orbifold Euler characteristic is defined as 

,».(a, = g(-i,'i;^. 

Or, if Q is the set of all the graphs Gki and v{G) is the number of vertices of the 
graph G, then we can write 



(5.2) x"'(C*) = J2 



ceo 
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Let ci"' be the graph complex where C^"' is spanned by aU graphs of valence 3 
or more on k vertices such that the number of edges minus the number of vertices 
is n — 1. (If the graph is connected, this means that there are n independent loops 
in the graph, i.e., the fundamental group of the graph has rank n.) We want to 
compute the generating function 

(5.3) x°'^'(a)(t) = ^x"^d"^)r 

n>2 

of the orbifold Euler characteristic of this graph complex. 
First, consider the power series 

(5.4) h{t,x,y) = eM-\E^^ixy)), 

and replace each occurrence of y'' by the number Chk (note that this eliminates 
all the odd powers of x). We get a power series in and x, which is almost the 
same as the series |T^|(a;) = (C/i x Par^^){x) considered in Section 3, but with an 
extra —t~^ for each vertex of the graphs. More precisely, the coefficient of x'^'^t~^ 
is the sum (—1)" J2 l/l-'^ut(G')| over graphs with e edges and v vertices. 
Now, substitute x = t^^^ to obtain the power series 

(5.5) h{t,y)=exp{--^E^+{t'^'y)). 

If wc siibstitutc Chk for each occurrence of y'^, we get a power series in t, where the 
coefficient of is the sum ^(— l)^("-^)/|Aut(G)| over graphs with n = #{edges} — 
#{vertices}. Note that = Ef+ = = ensures that all the exponents of t 
arc non-negative, and that if wc order h(t, y) by exponents of t, then the coefficient 
of is a polynomial in y, with exponents between 2n and 6(n — 1). 

So we have proved the following: if we replace each occurrence of y'^ by Ch^ 
in the power series h{t,y) defined above in (5.5), we get times the generating 
function X°'''(C*)(i)- 

Now it is a simple matter to check that the conditions of the Corollary to The- 
orem 1 are satisfied, and hence 



(5.6) 



1 /'OO 

-= / e-y"^\{t,y)dy^t-\°'^\C,){t) ast^O. 



The same method works for many different types of graph complexes. 

Theorem 2. Consider the graph complex C, consisting of graphs with some struc- 
ture at each vertex, represented by the species Q, where Ck is spanned by all such 
Q-graphs on k vertices. If Qo = Qi = Q2 = 0, then the orbifold Euler character- 
istic of the graph complex can be obtained from the power series 

(5.7) h{t,y)=eM-]Q{t'^^y)) 

by replacing each occurrence of y^ by Chk . 

//, moreover, there is a function (also denoted by h{t,y)) which has the asymp- 
totic expansion ( 5. 7) and which satisfies the conditions of the Corollary to Theo- 
rem 1, then the orbifold Euler characteristic is an asymptotic expansion: 

(5.8) -= / e-y''^h{t, y) dy « t-ix"^a)(i) ast^Q. 
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6. The three graph homologies 

In [10], Kontsevich constructed three graph complexes, spanned by what we wiU 
call commutative, associative and Lie graphs. The homology of the chain complex, 
i.e., the quotient ker(9)/im(9) is Kontsevich's graph homology. 

Commutative graphs are the simplest: they are finite graphs of valence 3 or 
higher, i.e., of the type considered so far. The species Q at each vertex is E^~^. 
Classes in commutative graph homology play a role in the theory of finite type 
invariants of homology 3-spheres [2] . 

Associative graphs have the additional structure of a cyclic ordering of the half- 
edges at each vertex. These are the same as the ribbon graphs or fat graphs of 
Teichmiiller theory. Associative graph homology is the direct sum of the homology 
of certain moduli spaces of mapping class groups [5] . The species at each vertex is 
C, the species of cyclic ordcrings; its terms are Cn — {n — 1)!. 

Lie graphs have a planar trivalent tree at each vertex, subject to some relations 
corresponding to the anti-symmetry and Jacobi identity of Lie algebras. The cor- 
responding graph homology compiitcs the homology (with trivial real coefficients) 
of the group Out(_F,i). The details can be found in [5]; for now we only need the 
number (more precisely: the dimension of the quotient vector space) of such trees 
with n leaves, which is LiCn — (n — 2)!. 

Other types of graph homologies can also be considered, and Kontsevich's ma- 
chinery, including the Euler characteristic computations explained here, can be 
applied to them; but only these three are known so far to have applications in 
other parts of mathematics. 

There are two complications. First, Kontsevich considers connected graphs only, 
thereby obtaining the primitive homology of the corresponding topological object. 
This is easy to remedy: the logarithm of the generating function for the orbifold 
Euler characteristic of all graphs yields that of connected graphs only. 

The second problem is that Kontsevich's graphs come with a sense of orientation, 
and the chain groups are spanned by those graphs only which do not have any 
orientation-reversing automorphisms. By a lucky coincidence, in the associative 
and Lie cases, removing the graphs with orientation-reversing automorphisms does 
not change the homology or the orbifold Euler characteristic; this is shown in 
[5]. In the commutative case, unfortunately, it does: already in the case of 3 
independent loops, we get x°''^ = if we include graphs with orientation-reversing 
automorphisms, and x""^^ = —1/48 if we do not. According to Kontsevich, we do 
get an approximately correct answer, though, because the typical graph does not 
have any automorphisms. 

Table 1 shows the orbifold Euler characteristics Xc'^^, Xa'^^ s.nd xf"^ of Kont- 
sevich's commutative, associative and Lie graph complexes, respectively, for small 
numbers of independent loops. The results are after taking the logarithm, but with- 
out removing graphs with orientation-reversing automorphisms. It was computed 
using (5.7) and the computer algebra software Mathematica. 

A striking feature of this "experimental data" is the equality of the commutative 
and the associative orbifold Euler characteristics. Furthermore, if we type these 
numbers into the On-Line Encyclopedia of Integer Sequences [14], we find that 

(6.1) (Xr')n = (Xr)n = -J^y 

n[n — 1) 

where the are the Bernoulli numbers. By using asymptotic expansions, we can 
prove these facts. 

Remark. The numbers in Table 1 are not new. The associative case (including (6.1) 
for Xo'^'') follows from Harer and Zagier's computations [8] of the Euler characteristic 
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Table 1 . The orbifold Euler characteristic of Kontsevich's graph complexes 





orbifc 


)ld Euler characteris 


3tic 


# loops 


Xc (commutative) 


Xa (associative) 


Xi (Lie) 




1 


1 


1 


2 


12 


12 


~24 


3 








1 

~48 




1 


1 


161 


4 


~360 


~360 


5760 


5 








367 

5760 


6 


1 


1 


120257 


1260 


1260 


580608 


7 








39793 
45360 


8 


1 


1 


6389072441 




1680 


1680 


1393459200 


9 








993607187 

34836480 


10 


1 

1188 


1 

1188 


5048071877071 
24524881920 


11 








9718190078959 
5748019200 



of mapping class groups; the power series xT'"^^) computed by Smillie and 
Vogtmann [16] as the generating function of the orbifold Euler characteristic of 
Out(_F,i). Our paper confirms their computations by describing another, simpler, 
derivation of these results. 



7. Evaluating the integrals 

The following computation in the commutative case is the same as Penner's [15] . 
The other two sections are (unfortunately incomplete) attempts at extending his 
method to the associative and Lie graph complexes, as well. 

7.1. The commutative case. The commutative series Q{x) = E^^{x) is the 
Taylor expansion of f{x) — ~ 1 — x — /2 around 0. We can apply Theorem 2 
directly with h{t,y) = exp(— i/(t^/^y)), which gives us 

(7.1) 

i"'x°'^0«log(^ P e-^'/2gxp(^-l(e^'^-l~2/Vt-2/V2)) dy^ 



oo 



= log ( ^ I e ete* dy 



/2'JTt 
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Here we substituted u = y\ft. Now substitute z = e"/i to obtain 

where in the last step we apphed the integral formula r(a;) = e~'^z^~^ dz which 
is one of several equivalent definitions of the gamma function. It is well known that 
r{z) has the asymptotic expansion around infinity given by Stirling's formula 

(7.3) r(^) « y^(2/e)^e^(^) as z ^ oo, 

(see, for example, [1]), where 

oo J-. / 1 \ 2n— 1 

and the arc the Bernoulli numbers. If we apply this formula to z = 1/t, take 
the logarithm, and compare terms, we get 

n—l ^ ' 

In other words, the orbifold Euler characteristic of ci™^ is zero for m odd, and it 
is given by (7.5) for m = 2n. Note that our notation for the Bernoulli numbers 

assumes 

(7.6) {Bi, ^2, Bs, -B4, -B5, . . .} = {-i i, 0, 0, . . .}. 
This seems to be standard, but it is different from the notation in [1]. 

7.2. The associative case. The associative series Q{x) = C^'^{x) is the Taylor 
expansion of f{x) = — log(l— x)— a;— around 0, which is undefined when a; > 1. 
However, the Taylor expansion of f(ix) is Q{ix), and h{t,y) = exp(— j/(it^/^t/)) 
satisfies the conditions of the Corollary to Theorem 1. Since the odd powers of 
it^/'^y are eliminated. Theorem 2 yields 

(7.7) log (-1= y ^ e-^V2ei/(»v7.) ^ r^'^^-t). 
Now we can write 



(7.8) 



= log (^-^ p (1 - iyVt)-"-e-'yl^'dy^ , 



where the logarithm function is defined on the complex plane minus the positive 
real axis. Next, we substitute z — j{iyVi— 1) to get 

i-'xr(-i) « / - (-te)-te-^-t dz 



(7.9) 

= log(^(et)-^ J^{-z)-h-^dz 
where 7 is the curve which comes from +00 above the positive real axis, goes 
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Y 



Figure 2. The contour 7 



around the origin counter-clockwise, and goes back to +00 below the real axis. 
This contour integral is related to the gamma function via Hankel's formula 

(7.10) / (-z)"~^e"^ dz = -2i sm{TTu)T{u); 

see, for example, [17], section 12.22. Applying this formula, together with the 
identity 

(7.11) ^T{u)V{l-u) 



sm ITU 



to our integral, we get 
(7.12) 



t-^-^-t) « log ( ^(e<)-^ sin Ui-^\\t(i-}- 



log ( V2TTt{ety- ^ 



r(iA). 

This is exactly the negative of the expression (7.2) we got in the commutative 
case! Hence, if we substitute —t back in place of t, we get the same orbifold Euler 
characteristic 



2n 



00 

n— 1 ^ ^ 

for the associative graph complex as for the commutative graph complex. (Which, 
however, is not the exact answer in the commutative case, as we have noted earlier.) 

7.3. The Lie case. I do not know yet how to deal with this case. The Lie series 
Q{x) = Lie{x) is the Taylor expansion of f{x) — {1 — x) log(l — x) — 1 + x — x^/2 
around 0, which is once again not defined for x > 1, but f{ix) is. Therefore, in the 
same way as in the associative case, from Theorem 2 we get 

(7.14) 



/2tt J- 



l{{l-iy^)log{l-iyVi)-l+iyVi) 



= log ( ^ / g..log(^)-^+^„log(^) ] ^ 



where w = j{l — iy^/i). If we now substitute 
(7.15) t-\r\~t) « log 



SR(z) 




I do not know whether this integral can be related to some known asymptotic 
series, like in the other two cases; but even if it can not, we might be able to use it 
to extract some information about the coefficients of x"*^ (t) ■ 
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8. Summary and further questions 

We have derived Kontsevich's formula (5.8) using the method of Feynman dia- 
grams and the theory of asymptotic expansions. Using these techniques, one could 
also compute the orbifold Euler characteristic of other graph complexes, as stated 
in Theorem 2. 

I have written a Mathematica notebook which performs this computation; it can 
be downloaded from http://www.math.cornell.edu/Research/Dissertations/ 
Gerlits/code . One enters the formula for on the second line of the note- 

book, executes each of the subsequent lines, and the output of the last line is the 
orbifold Euler characteristic of the Q-graph complex for 2 < n < 11 independent 
loops. The numbers in Table 1 were generated by this notebook. 

For example, one could consider the "chord diagram graph complex" by putting 
a chord diagram at each vertex. Then Q = Ch^~^, and the generating function of 
the orbifold Euler characteristic is 

So To 131 4 449 , 80179 . 16459 , 127239605 » 

- H H f H r 

8 16 128 128 5120 192 229376 

16956565 g 27521691751 6769184257 „ 
4096 786432 20480 

Kontscvic;h claims that in the Lie case, one gets even larger numbers for the orb- 
ifold Euler characteristic than the Bernoulli numbers which show up in the other 
two cases. This ought to be able to be proved using (7.15); also, one should be able 
to give a direct proof for the fact that the orbifold Euler characteristic of the asso- 
ciative graph complex is zero when the number of loops is odd. Another direction 
would be to try to incorporate Kontsevich's orientation into the Feynman diagram 
apparatus so that we could exclude graphs with orientation-reversing automor- 
phisms, and compute the actual orbifold Euler characteristic in the commutative 
CcLSG) cLS well as for graph complexes based on other species. 
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